Violations of Bell's inequality for Gaussian states with homodyne detection and nonlinear 

interactions 
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We show that homodyne measurements can be used to demonstrate violations of Bell's inequality with Gaus- 
sian states, when the local rotations used for these types of tests are implemented using nonlinear unitary op- 
erations. We reveal that the local structure of the Gaussian state under scrutiny is crucial in the performance 
of the test. The effects of finite detection efficiency is thoroughly studied and shown to only mildly affect the 
revelation of Bell violations. We speculate that our approach may be extended to other applications such as 
entanglement distillation where local operations are necessary elements besides quantum entanglement. 



Violation of Bell's inequality |jj,[2[], which means failure of 
local realism, is perhaps the most profound yet controversial 
feature of quantum mechanics. It was Einstein, Podolsky and 
Rosen (EPR)'s work which challenged the completeness of 
quantum mechanics as a theory [3]. Although the original 
formulation of EPR's paradox involved the state of a bipartite 
system having a continuous spectrum, Bohm's version of the 
problem [4] and the seminal work by Bell [1, 2] moved the 
debate towards its discrete version, which quickly became the 
paradigm in the physics community. 

The proved experimental handiness of continuous variable 
(CV) systems, epitomized by proposals and realizations of 
schemes for quantum teleportation [5], among other exam- 
ples, has redirected considerable interest towards the investi- 
gation of Bell's inequality with these states. It has been known 
since Bell that the original EPR state allows a realistic de- 
scription in terms of the canonically conjugated variables of 
position and momentum because its Gaussian Wigner func- 
tion can be used as a classical probability distribution for a 
hidden variable model [2]. Later, however, it has been proven 
that a two-mode squeezed vacuum state, whose Wigner func- 
tion is a Gaussian, can violate Bell's inequality, although non- 
optimal, in the joint measurement of phase-space displaced 
parity operator [6]. The loophole here is that the inclusion 
of photon counting measurements (necessary for parity de- 
termination) negates a reakistic interpretation of the Wigner 
function [7]. This result has triggered extensive investigation 
on the interplay between measurements and Gaussian charac- 
ter of CV states in Bell's inequality tests. It has been found 
that non-Gaussian CV entangled states can be used to demon- 
strate violations of Bell's inequality by means of Gaussian 
measurements (i.e. measurements that preserve the Gaussian 
nature of a state, such as homodyne detection) [8]. It has also 
been shown that a non-deterministic Bell's inequality test can 
be devised using Gaussian CV states to show Bell violations 
when non-Gaussian conditioning measurements are combined 
with homodyne measurements |2]. However, in virtue of 
Bell's argument for a realistic description of an EPR state, it 
is well known that Bell's inequality tests in which both states 
and measurements are Gaussian, are destined to satisfaction 
of Bell's inequality [2]. More practically, this "no-go" result 



seems paired with the impossibility of entanglement distil- 
lation for Gaussian states under Gaussian-preserving opera- 
tions II10I1 . which in turn limits the implementation of efficient 
quantum repeaters for Gaussian states 11X111 . 

In this paper, we describe a different angle on the afore- 
mentioned problems. We examine Clauser-Horne-Shimony- 
Holt (CHSH)'s version of Bell's inequality 11211 and show its 
violation by Gaussian states subjected to nonlinear local uni- 
tary operations and homodyne measurements. In order to effi- 
ciently illustrate our findings, which we describe with respect 
to routinely generated entangled Gaussian states, it is conve- 
nient to first address a situation where the local operations re- 
quired for our Bell's inequality test are treated as ideal single 
qubit rotations [13] for coherent-state qubits [ 14, 15]. We then 
present physical transformations IU4ll which are able to repro- 
duce the core features of the idealized case (cfr. Fig. [TJ. In 
both ideal and physical cases, a considerable (although non- 
optimal) degree of violation of the Bell-CHSH test is found, 
resilient to imperfections at the detection stage. Our study 
shows that the local unitary operations, the importance of 
which have not previously been carefully examined, may play 
a key role in utilizing Gaussian states for Bell's inequality tests 
and, possibly, other applications such as entanglement distil- 
lation. Our protocol is shown to be quite resilient to spoiling 
effects such as non-zero temperatures and low detection effi- 
ciency. We also stress that our scheme makes use of Gaus- 
sian squeezed states, which are resources routinely employed 
in all-optical experiments dealing with continuous variables. 
This distinguishes our proposal from the case of entangled co- 
herent states being used for Bell's inequality tests (such as in 
Ref. fl3IO . which are non-Gaussian and more demanding to 
produce. 

The remainder of this paper is organized as follows. In Sec- 
tion|T]we describe the formal approach to our Bell's inequality 
test. We introduce the Gaussian states we probe and the class 
of nonlinear unitary operations we consider. This is achieved 
by first studying a sort of idealized case and then moving to 
the real physical situation. We show that Bell's inequality can 
be violated by Gaussian states with homodyne measurements 
using the specific class of operations. In Section HT1 we study 
the effects of detection inefficiencies and show that they can 



be counteracted by increasing the squeezing of the initial re- 
source. This very same strategy can be used in order to cope 
with a mixedness initial state instead of a pure, ideal resource. 
This is consistent with a previous result 1 16] using entangled 
thermal states [17], where the same effects could be achieved 
by increasing the distance between the component thermal 
states [16]. Section|lII]shows that the required level of squeez- 
ing to show Bell violations can be considerably reduced when 
another class of experimentally relevant Gaussian states are 
used. Finally, in Section [TV] we summarize our results and 
discuss their physical implications. 



I. TEST FOR BELL'S INEQUALITY 



superposition of coherent states 11 



In Ref. 111311 . it was shown that a superposition of two co- 
herent states, |c_|_) = N+(\a) + \ — a}) with the normalization 
factor N + and coherent states | ± a) of amplitudes ±a, when 
divided at a beam splitter, violates Bell-CHSH inequality us- 
ing homodyne measurements and nonlinear interactions. One 
can show that the fidelity T between a coherent-state super- 
position |c+) and a single-mode squeezed state is very high 
when a is relatively small (e.g. T > 0.99 for a < 0.75). This 
motivates us to first investigate violation of Bell's inequality 
for single-mode squeezed states divided at a beam splitter. We 
shall later study another set of Gaussian states which outper- 
form the results for this case. 

Let us suppose that two parties, Alice and Bob, share an en- 
tangled state generated using a single-mode squeezed vacuum 
and a 50 : 50 beam splitter [18]. Analytically, the state can 
be described by the following Gaussian-weighted continuous 
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FIG. 1: Schematic of a Bell-CHSH inequality test with Gaussian 
states and homodyne measurements. A quantum correlated two- 
mode Gaussian states is produced at a source and distributed to Alice 
and Bob who perform local effective rotations and homodyne mea- 
surements. The effective rotations are physically implemented by 
cascading Kerr-type nonlinearities and phase-space displacement op- 
erations, as shown in the inset. The Gaussian state produced by the 
source can be either a two-mode squeezed state or the state resulting 
from the superposition, at a 50 : 50 beam splitter, of vacuum and a 
single-mode squeezed state. 
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where Q(r,a) — exp[— (1 — tanhr)a 2 /(2tanhr)], r is the 
squeezing parameter, a £ M and J\f = 1/V27rsinhr is 
the normalization factor. The class of nonlinear transforma- 
tions we consider can be understood as an approximation of 
the following rotations performed in the bidimensional space 
spanned by the generic coherent state {|±/J}} (/? S C) 111311 
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where 6j is the effective "angle" of such idealized rotations 
and j = A, B labels Alice's or Bob's site. It should be noted 
that the "idealized" transformation described in Eq. (f5]l is not 
unitary (approximately unitary when (3 is large) so that it can- 
not be performed deterministically. The actual physical local 
transformation using nonlinear interactions will be considered 
later in this Section. 

After the application of the local operations (fJJ) to their 
respective mode, Alice and Bob perform bilocal homodyne 
measurements, which result in the joint probability-amplitude 
function 
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with | a:) (\y)) the in-phase quadrature eigenstate of Alice's 
(Bob's) mode. A sketch of the thought-experiment is pre- 
sented in Fig. 03 In order to test CHSH version of Bell's in- 
equality, we need to construct a set of bounded dichotomic ob- 
servables, which we do by assigning value +1 to a homodyne - 
measurement's outcome larger than 0, and —1 otherwise [8]. 
With this, a joint probability of outcomes can be calculated as 
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(4) 



where the subscripts fc, I = ± correspond to Alice's and 
Bob's assigned measurements outcomes ±1 and the integra- 
tion limits are such that + s = oo, +/ = — s = and 
— i = — oo. We can now calculate the Bell-CHSH function, 
B(e A ,0 B ,6' A , ff B ) = C(0 A , B ) + C($' A , B ) + C(0 A , o' B ) - 
C(9' A , 9' B ), where we have introduced the correlation function 
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According to local-realistic theories, the Bell-CHSH inequal- 
ity \B(6 A ,0 B ,0' A ,0' B )\ < 2 holds. Quantitatively, we have 
found that 

2arctan(sinhr)cos(40,4)cos(40.B) ,„ 

(->id(0 A ,0B,'r') = . ,.- , (6) 

1 7 tt(1+ E sin(40 i )[22^+sinhr-]) 

with j, k = A,B and the subscript id is used in order to re- 
mind of the idealized version of local operations being used. 
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FIG. 2: Bell test for ideal rotations. The Bell-CHSH function is plot- 
ted against the squeezing parameter r for three values of the detection 
efficiency r\. We show the case corresponding to ideal homodyne de- 
tection (solid line), r\ = 0.5 (dashed line) and r\ = 0.05 (dotted line). 
The horizontal line shows the bound for local realistic theories. 



The behavior of the numerically optimized Bell-CHSH func- 
tion corresponding to Eq. © is shown by the solid curve in 
Fig. |2 which demonstrates that a local realistic description of 
10 ab * s impossible as the squeezing parameter for the initial 
single-mode state surpasses ~ 2.1. The degree of violation of 
the Bell-CHSH inequality then reaches a maximum of ~ 2.23, 
robustly against r. 

Now that we have gained a quantitative picture of the be- 
havior of the Bell function under the class of formal opera- 
tions and Gaussian measurements considered in our work, it 
is time to provide a physically effective description of each 
rotation Rj(9j). Such physical implementation stems from 
the observation made in Ref. [13] that Eqs. (O can be ap- 
proximated by a combination of single-mode Kerr interac- 
tion U Ken = e _i&KOT * with H Ken = Ml(tfa) 2 (Q being the 
strength of the non-linear coupling) and displacement of am- 
plitude <p G C, D(<p) = e v^-v*&_ Here a (a*) i s the an- 
nihilation (creation) operator of a field mode. A single-mode 
Kerr interaction may be implemented, for example, by nonlin- 
ear crystals [20, 21] while the displacement can be easily per- 
formed via a beam splitter with high transmittivity and a local 
oscillator. In detail, the evolution induced by the effective ro- 
tations Vj(9j) = L r KeiT-D(*^i/rf)C^Ken- on input coherent states 
|±/3) ■ is given by the following expressions ((3 = j3 r 
and dgt determines the amplitde of the displacement) [ 

V ] {0 ] )\P)^\[e^^{\P+ l -^) ] + l \-P- l -^) ] ) 

+fe-4/»r(|-/3+^> J +<|/J-!^> J )} 1 
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tary operation. The physical operation Vj(6j) is a good ap- 
proximation of the ideal operation Rj(9j) when the ampli- 
tudes of coherent states on which the operation is acted are 
large. As seen in Eq. (fl}, our squeezed state can be expanded 
in terms of coherent states with a Gaussian weight factor as 
a function of the coherent amplitude. When the squeezing r 
is large, contributions of coherent states of small amplitudes 
will become arbitrarily small. This implies that as the squeez- 
ing r becomes large, the results of the Bell-CHSH inequality 
violation using the ideal rotation Rj (9j ) should be closer to 
the results using the physical rotation Vj (9j). 

We then adjust our notation and indicate with 

C ef {9 A ,9 B ,x,y) = \(x,y\V A (9 A )V B {9 B )\0 AB \ 2 the 
probability of measuring the values x and y of the quadra- 
ture variables at the homodyne detectors. The subscript 
clearly states that this is the function associated with the 
use of physical effective rotations. Quantitatively, C e f 
is easily found using the projection of a coherent state 
onto a position quadrature eigenstate \x), which is given 
by (x\P) = Tr-i/^v^ftz-i^-v^,) 2 -^, (18 ^ We 
eventually obtain 
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Note that Vj(9j) is unitary while Rj(9j) is not strictly a uni- 



(8) 

We are now in a position to build up the Bell-CHSH func- 
tion for our Bell's inequality test in such physically effective 
case. Unfortunately, producing an analytic result is rather de- 
manding due to the semi-infinite range of integrations over 
the quadrature variables x and y, which also enter into the 
trigonometric functions in Eq. (|S), required in order to gather 
the joint probabilities Pki(9 A ,9 B ). We have therefore per- 
formed the Bell's inequality test by numerically evaluating 
the Bell-CHSH function for a set value of d and by scan- 
ning the squeezing parameter r. The results are shown by 
the top-most curve in Fig. [3] where violation of local realis- 
tic theories starting from r > 2.1 is observed, which is in 
full agreement with the ideal-rotation case. Also, the degree 
of violation is consistent between the two cases, |-B| maa; be- 
ing 2.229 at r = 3.3. Although the reproduction of the be- 
haviour for large r is computationally demanding, it is pos- 
sible to perform a qualitative comparison between ideal and 
effective case by looking at the corresponding joint probabil- 
ity functions \Cid{9 A , 9 B ,x, y)\ 2 and Eq. ([8]), evaluated at the 
angles corresponding to the (numerically-optimized) associ- 
ated Bell-CHSH function. This is done in Fig. |4] where the 
clear similarity of the two probability functions ensures the 
closeness of the value of the corresponding Bell-CHSH func- 
tions. 



II. ROBUSTNESS TO IMPERFECTIONS 

Although homodyne detectors have rather high efficiencies, 
the violation of the Bell-CHSH inequality by \£,) AB i s far from 



2\/2, the maximum given by Tsirelson's bound [22]. One 
might thus wonder whether even mild detection inefficiencies 
are sufficient to wash out the Bell-CHSH inequality violation 
unveiled in Fig. [3] An important issue to address is thus given 
by the effects of detection inefficiencies. As done before, we 
first gain an idea of the expected behavior by studying the 
idealized picture. 

In order to quantitatively assess this point, we have mod- 
eled the imperfect homodyne detector onto which mode j = 
A, B impinges as the cascade of a beam splitter of transmit- 
tivity 77, mixing mode j to an ancillary vacuum mode aj, 
and a perfect homodyner. We are not interested in the state 
of the ancillae, which are discarded by tracing them out of 
the overall state, so that \Cid(9 A , &B,x,y)\ 2 is changed into 
(x,y\Tr aAaB ip(9 A ,6 B )\x,y) with 




FIG. 4: We compare the behavior of the joint-probability functions 
\Cid{O A ,0B,x,y)\ 2 and \C e f(dA,6 B ,x,y)\ 2 against the quadrature 
variables x and y for r — 4. The angles 8a, b are those maximizing 
the corresponding Bell-CHSH function. We thus have 9a = 0.061 
and 6 B = 0.182 (6a = -0.009 and 6 B = 0.004) for the left- 
most (rightmost) plot. Moreover, f J dxdy\C e f (6a, 6b, x,y)\ 2 ~ 
J f dxdy\dd(6A, 6b, x, y)\ 2 , regardless of the domain of integra- 
tion. 
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where the beam splitter operation between modes j and 
the corresponding ancilla aj is defined as Bj a (£) = 

exp[|(a]& aj — cijbl)] with cosC = ^/r\ and b aj being the 
annihilation operator of aj 112311 . The remaining procedure for 
the construction of the appropriate Bell-CHSH function re- 
mains as described above. The final form of the correlation 
function, which now depends on the efficiency as well, is ob- 
tained from Eq. (|6]l by simply replacing arctan(sinhr) — * 
arctanf /1 T |'l M " r , ). The behavior of the associated Bell 

x V l-\-2r/e r sinn r y 

function is shown, for two values of rj, in Fig. [2] We ob- 
serve a rather striking robustness of the Bell function with re- 
spect to the homodyners' inefficiency: Even severely ineffi- 
cient homodyne detectors would be able to unveil Bell-CHSH 
inequality violations with a state which is initially squeezed 
enough. By simply increasing the squeezing of the input state, 
one can compensate the effects of detection inefficiencies. Al- 
though for small values of 77, the required squeezing factor be- 




comes prohibitively large, the trend revealed by the ideal case 
leaves quite a few hopes for the physical effective one as well. 
In fact, such robustness persists when the local operations (0 
are used, as shown in Fig. [3] for 77 = 0.8 and 0.3 (chosen for 
easiness of representation). The squeezing threshold at which 
the Bell's inequality test starts to be violated increases only 
quite slowly as the quality of the homodyne detectors is de- 
graded. In passing, we should stress that the beam-splitter 
model used for the description of an inefficient homodyne de- 
tector can be used in order to describe the influences of ex- 
ternal zero-temperature reservoirs coupled to the correlated 
two-mode state we are studying. Thus, similar conclusions 
regarding the resilience of the Bell-CHSH function to losses 
induced by a low-temperature environment can be drawn. 

We complete our study about the effects of imperfections 
by investigating the case in which we start with a mixed re- 
source. This is practically quite relevant, given the fact that, 
experimentally, a single-mode squeezed thermal state is in 
general produced instead of a pure single-mode squeezed vac- 
uum state. This is formally accounted for by considering the 
resource state 



P 8 1b- 
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a a a a 

(10) 
where Tin, a) = e~l Q ~ d l ' n /im. (a — a r + laA is the 
Glauber-Sudarshan function of a single-mode state at ther- 
mal equilibrium with mean photon number ri and displaced, 
in phase space, by d G R JUt] while S A (r) = e^ (a 2 -°- 2 ) 
is mode-A squeezing operator. Eq. ( TTOb results from super- 
imposing at a 50 : 50 beam splitter a squeezed displaced 
thermal state of mode A and the vacuum state of mode B. 
It is straightforward to find that Eq. ( TTOb can be written 
™Pab = Jd 2 af{r,V,a)\a/V2,a/V2) AB (a/V2,a/V2\ 
with V = 2n+l and 



FIG. 3: The numerically optimized Bell function is plotted against 
the squeezing parameter r for the case of physical effective rotations 
and three values of detection efficiency. The horizontal line shows 
the bound for local realistic theories. The actual value of d is ir- 
relevant, in this figure. The solid line with rj = 1 embodies the 
ideal-detector case with the other two cases of n = 0.8 and rj = 0.3. 
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This state is then locally rotated and projected onto quadra- 
ture eigenstates by means of homodyne measurements. Once 
more, for clarity of our arguments, we refer to the case of ideal 




FIG. 5: Bell-CHSH test for an input squeezed thermal state superim- 
posed to vacuum at a 50 : 50 beam splitter. The Bell-CHSH function 
is plotted against r and V = 2n + 1, i.e. the thermal variance of 
the state. The horizontal plane shows the bound for local realistic 
theories. 



rotations. The use of our formal procedure applied so far lead 
to the correlation function 
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(12) 
Clearly, C s t = Cid when V = 1, i.e. when a pure state is 
generated. With this expression, one can easily build up the 
Bell-CHSH function and test its behavior against the thermal 
parameter V and, as usual, the squeezing. The results are 
shown in Fig. |5J where it is shown that it is enough to consider 
a slightly more squeezed initial resource in order to counter- 
act any thermal effect. The same conclusions are reached by 
using the set of nonlinear unitary transformations Vj(9j), al- 
though the analysis is largely numerical and more involved. 



smaller values of r are taken. Our starting point is the obser- 
vation Sill 



B 
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where we have used the single-mode squeezing operator 
Sj(r) — exp[^(a3 — a"- )] (j = A, B) and its two-mode 
version S AB {r) = exp[r(a A a B — a A a B )]. Therefore, our re- 
source \0ab i s f° rma lly equivalent to a two-mode squeezed 
state that is also subjected to additional local squeezing oper- 
ation. These latter are unable to change the nonlocal content 
of the state being used and could well be regarded as a pre- 
stage of the local actions (comprising nonlinear rotations and 
homodyne measurements) performed at Alice's and Bob's site 
respectively. We now remove them from the overall setup for 
Bell's inequality tests by considering, instead of Eq. (Q]i, the 
standard two-mode squeezed vacuum [26J 



\?)=M j 'd 2 f3G'(r,l3)\P,f3* 



with weight function (|27| 



g'(r,0)=exp[ 



1 — tanhr| 
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(15) 



and normalization factor M. — (tv sinhr) -1 . The adaptation 
of the formal procedure described in our work to the use of 
this Gaussian resource is quite straightforward. For the sim- 
ple case of ideal local rotations, the correlation function for 
joint outcomes at Alice's and Bob's site is identical to Eq. (O 
with the replacement r — * 2r. The violation of the local real- 
istic bound occurs now for r ~ 1 and the entire Bell-CHSH 
function shown in Figs. [2] (a) is "shifted back" on the r axis 
accordingly This effect is the same when Eqs. (0 are used, al- 
though the form of C e / is too cumbersome to be shown here. 
For clarity, we note that a two-mode squeezed state of squeez- 
ing r can be generated using two single-mode squeezed states 
of the same degree of squeezing and a beam splitter as 



III. IMPROVEMENT USING TWO-MODE SQUEEZED 
STATES 

The required level of squeezing, e.g. r > 2 for r/ > 0.8, re- 
vealed in Figs. 2 and 3 to demonstrate Bell-CHSH inequality 
violations is experimentally difficult to achieve using current 
technology. In this Section, we show that this requirement can 
be radically reduced by using another class of Gaussian states. 
So far, we have investigated the Bell's inequality test under 
nonlinear operations using the paradigmatic source given by 
state \£,)ab- However, the behaviour of a Bell-CHSH function 
strongly depends on intrinsic properties of the tested quan- 
tum correlated state. In fact, this can be seen as the "dual" 
of the well-known fact that the same bipartite entangled state 
behaves differently, in terms of Bell inequality tests, under 
different sets of local operations. Here, we are interested in 
finding out whether another realistic Gaussian resource is con- 
ceivable for the violation of our Bell-CHSH inequality when 



S AB (r) |00) AB = B AB (-)S A (r)S B (-r) |0,0) 
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This means that single mode squeezed states of \r\ > 1 
(> 8.7dB) can be used as resources to show violations of 
Bell's inequality. This makes our proposal closer to an ex- 
perimental implementation as such high levels of squeezing 
can be generated (for example, up to lOdB [28]) using cur- 
rent technology. On the other hand, the local nonlinear oper- 
ations may be more demanding and various types of unitary 
interactions need to be investigated to improve experimental 
feasibility of our approach. 



IV. CONCLUSIONS 

We have shown a way to unveil violations of Bell's inequal- 
ity for two-mode Gaussian states by means of nonlinear local 
operations and Gaussian homodyne measurements. Besides 



its theoretical interest, which stays at the center of current in- 
vestigations on entangled CV systems and their fundamental 
features, our study emerges as an appealing alternative to the 
current strategy for Bell's inequality tests based on the use of 
appropriately de-Gaussified resources and high-efficiency ho- 
modyning. Our proposal has been shown to be robust against 
the inefficiency of the homodyne detection and mixedness in 
the initial resource. This robustness is consistent with a previ- 
ous study using entangled thermal states [16]. 

While the squeezing degree of r > 1 required for resource 
Gaussian states is possible to achieve using present day tech- 
nology, the strong nonlinear interactions required to imple- 
ment the local operations may be more difficult to realize. On 
the other hand, it is worth noting that there has been remark- 
able progress to obtain strong nonlinear effects [21, 29]. 

There remains some interesting future work. As the local 
operations used in our paper are not necessarily optimized 
ones, the research for more efficient local operations is de- 
sirable. Since we have used nonlinear operations to reveal vi- 



olations of Bell's inequality for Gaussian states and Gaussian 
measurements, it is also natural to extend this investigation to 
entanglement distillation protocols for the Gaussian states. In- 
teresting open questions are therefore whether there exist such 
entanglement distillation protocols using the type of local op- 
erations employed in this paper and how much they would be 
feasible and useful. 
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